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We present the first cosmological measurement derived from a galaxy density field subject to 
a ‘clipping’ transformation. By enforcing an upper bound on the galaxy number density held in 
the Galaxy and Mass Assembly survey (GAMA), contributions from the nonlinear processes of 
virialisation and galaxy bias are greatly reduced. This leads to a galaxy power spectrum which is 
easier to model, without calibration from numerical simulations. 

We develop a theoretical model for the power spectrum of a clipped held in redshift space, which 
is exact for the case of anisotropic Gaussian helds. Glipping is found to extend the applicability 
of the conventional Kaiser prescription by more than a factor of three in wavenumber, or a factor 
of thirty in terms of the number of Fourier modes. By modelling the galaxy power spectrum on 
scales k < 0.3/iMpc“^ and density huctuations < 4 we measure the normalised growth rate 
fasiz = 0.18) = 0.29 ±0.10. 


I. INTRODUCTION 

The spatial distribution of galaxies encodes a wealth of 
information relating to the composition and evolution of 
the Universe. The apparent positions of galaxies in red¬ 
shift space offers a glimpse into both the density and ve¬ 
locity perturbations associated with dark matter. These 
in turn are inhuenced by a number of phenomena in fun¬ 
damental physics, such as the mass of the neutrino and 
the nature of gravity. However two key factors have thus 
far restricted our view: (a) the advanced stages of grav¬ 
itational collapse are highly unpredictable and (b) the 
uncertainty associated with galaxy bias, defined as the 
manner in which the galaxy distribution reflects the dark 
matter distribution. In Fourier space, conventional anal¬ 
yses impose a maximum wavenumber beyond which the 
data points are considered unpredictable and are sim¬ 
ply discarded. For example, despite utilising numerical 
simulations to calibrate the non-linear power spectrum, 
recent studies of redshift space distortions typically trun¬ 
cate the power spectrum at /c < 0.2 hMpc~^ [1] or exclude 
galaxy pairs closer than 25/i“^Mpc [2]. While most of 
the nonlinear behaviour is successfully disposed of, so too 
is much of the cosmological information. 

A number of different methods have been proposed to 
allow smaller clustering scales to be exploited. A phe- 
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nomenological model has been developed by Kwan et al. 
[3], and expanded by Linder and Samsing [4], to model 
the power spectrum based on fits to numerical simula¬ 
tions. Reid et al. [5] extract information from small-scale 
clustering in BOSS using a model for the halo occupa¬ 
tion distribution. While achieving a significant increase 
in precision, this technique also relies heavily upon cali¬ 
bration from numerical simulations. Alternatively, vari¬ 
ous local transformations have been explored as a means 
of reducing the influence of nonlinearities, such as Gaus- 
sianisation [H [7] and the logarithmic transformation [8] . 
In Simpson et al. [9] it was shown that simply enforcing 
a maximum density could greatly increase the number 
of Fourier modes that could be modelled with the stan¬ 
dard set of theoretical tools. Recently it has also been 
found to enhance the observational signature in models 
of modified gravity which invoke a screening mechanism 
m- This clipping technique serves as the focus of the 
present work. 

Our inability to model the small scale power spectrum 
is not the only source of information loss. More funda¬ 
mentally, the onset of nonlinear gravitational collapse de¬ 
grades the amount of information held by the power spec¬ 
trum m-, this loss occurs for any non-Gaussian field. To 
extract some of this missing information, we can either 
perform additional measurements, such as higher-order 
statistics, or manipulate the field prior to evaluating the 
power spectrum. Previous attempts to extract cosmolog¬ 
ical information from the spatial distribution of galaxies 
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beyond the conventional two point statistics include mea¬ 
suring three-point statistics [HHni, Minkowksi function¬ 
als m, and the shapes of voids [16]. Local transforma¬ 
tions also appear promising in restoring this information 
to the power spectrum (see for example jHl ITTHT^ b Car- 
ron and Szapudi m present a transformation that is op¬ 
timised for extracting information from the power spec¬ 
trum of a Poisson-sampled field. Our approach differs 
slightly: instead of optimising the extraction of all infor¬ 
mation from the observable galaxy density field we 
seek to selectively extract only the predictable informa¬ 
tion from the field. At the highest values of pg both the 
nonlinear structure of the dark matter and the environ¬ 
mental impacts on galaxy formation present formidable 
obstacles in interpreting its value. Conversely, regions 
closer to the mean density are expected to behave in a 
more predictable and robust manner. The action of clip¬ 
ping preserves the location of galaxy clusters as useful 
information while discarding information relating to the 
precise value of their density contrast. 

Clipping is already known to be a highly effective 
technique for improving the theoretical modelling of the 
galaxy bispectrum [ 20 | and power spectrum [9| , when ap¬ 
plied to fields in real space. Before we can apply this tech¬ 
nique to data from real surveys, it must first be verified 
in redshift space. It would also be desirable to develop 
a deeper understanding of why it has been successful. 
These are two of the goals of this paper. The third is to 
apply clipping to the GAMA survey in order to obtain a 
low-redshift measurement of the normalised growth rate 
fcrg. 

In ^we review the theoretical background associated 
with the two point statistics of clipped fields and con¬ 
sider its extension to anisotropic fields. This theoretical 
framework is placed into a cosmological context in §III[ 
where we develop a model for the form of the clipped 
galaxy power spectrum. In ^TVjwe apply clipping to sim¬ 
ulated dark matter and galaxy fields in redshift space, 
with the results illustrated in ^V| The GAMA dataset 
is introduced in §VI[ w hile the main results of this work 
are presented in §VII| before our concluding remarks in 

^vini 

II. STATISTICAL PROPERTIES OF CLIPPED 
FIELDS 

Clipping is a local transformation characterised by the 
application of a saturation value (5o to a scalar field (5(x) 
such that 


and corresponding power spectrum Pc{k). We note that 
the clipping transformation induces a non-zero mean in 
Sc(x), which is why must be careful to specify the more 
general form of the autocorrelation function, as defined 
by Q. In the analysis of cosmological fields the sub¬ 
traction of the mean is conventionally omitted from this 
definition, since the mean is usually zero by construction. 

We begin by reviewing the special case where S(x) is an 
isotropic Gaussian field before generalising to anisotropic 
and higher-order fields. 


A. Isotropic Gaussian Fields 


For the case of Gaussian fields we may invoke Price’s the¬ 
orem [211 [ 22 ] to evaluate the two-point statistics associ¬ 
ated with the clipped field Sc(x) in terms of the original 
correlation function ^(r): 
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where is the local transformation defined by 0 ,?c{r) 
is the correlation function of the transformed field, and 
^2 5 '^) is the joint probability distribution for a Gaus¬ 
sian process 
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The functional derivative of the clipping transforma¬ 
tion is unity below the threshold, and zero above the 
threshold. This simple behaviour permits an analytic so¬ 
lution of the clipped correlation function ^ 0(^)5 which is 
given by |9] 


^c(r) = fUir) 


n=l 


C(0 


-I n+1 


Cn{uo), (5) 


where is the variance of the field prior to clipping, uq 
is the normalised threshold value uq = (^o/v^cr, is the 
fraction of the field which lies below the threshold 


/c = 2 [1 + (^ 0 )] > 

and Cn{x) is the distortion coefficient 

-^n—1(^) —2x^ 


Cn{x) = 


7T2^{n + 1)! 


( 6 ) 


(7) 


Jc(x) = So {S{x) > ^ 0 ) . . 

^c(x) = ^(x) (^(x) < So), 

yielding the clipped field Sc(x). In this section we ex¬ 
plore the statistical properties of Sc(x), with particular 
attention paid to its autocorrelation 

Cc(r) = ([(5c(x) - 4] [(5c(x + r) - 4]), (2) 


where Hn{x) is the Hermite polynomial of order n. Be¬ 
yond a scale-independent reduction in amplitude, clip¬ 
ping induces a distortion in the shape of the correlation 
function. However provided the slope of the spectral 
power is not too steep (|n| < 3), and the clipping remains 
weak {uo ^ 1 ), the leading order correction makes only a 
small contribution to the resultant power. Furthermore, 
terms at higher values of n decay rapidly. 
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B. Anisotropic Gaussian Fields 


For the more general case of anisotropic Gaussian fields 
we may express the two-point correlation function ^(r, 0) 
in terms of both the pair separation r and the orienta¬ 
tion vector 0. In deriving the expression given by ^ 
we made use of the joint probability distribution for a 
Gaussian process as given by (0, which is not directly 
applicable to anisotropic fields. However we can proceed 
by applying this single-parameter transformation sepa¬ 
rately at each fixed value of 6>, which in itself constitutes 
a single-parameter Gaussian process. Since the variance 
of each subspace is the same, = ^{r = 0,0) for all ori¬ 
entations, the transformation maintains the same func¬ 
tional form for all values of 0. Therefore the expression 
^ can be generalised to anisotropic fields. 

Transforming to Fourier space leaves us with the ex¬ 
pression for the clipped power spectrum 


Pe{k) = flP{k) + <72 ^ K) 

n=l 



>K(n-|-l) 


(8) 


where the notation represents a self-convolution of or¬ 
der n. In practice it is computationally more straightfor¬ 
ward to evaluate the higher order terms using powers of 
the correlation function, rather than performing multiple 
convolutions of the power spectrum. Further details of 
this calculation, as applied to mock cosmological density 
fields, can be found in Appendix [B] 

In line with the definition of the correlation function 
0, our definition of the power spectrum in ^ is speci¬ 
fied in terms of the mean subtracted field. For a clipped 
Gaussian field the mean is given by 

Sc = (5o(l - fc) - exp(-(5^/2(T2), ( 9 ) 

V ZTT 

In practice, it is not critical to account for this constant 
offset since it only contributes to the power spectrum at 
k = 0 . 


C. Second Order Anisotropic Fields 

In order to gain insight into how higher order terms 
respond to clipping, we repeat the procedure above using 


the square of a Gaussian random field, e{x) = S‘^{x) — 
<j^. The correlation function of this second order field, 
clipped at e = (^0 5 is well approximated by [9] 

Cc(r,e) = a2^(r,e)+e> ^ ’ 

(10) 

1 r '212 

a2 = - v^erf(i4o) - 2uQe , 

TT L 

(11) 

, V<5o + a-2 

(12) 


where a refers to the standard deviation associated with 
the original Gaussian field, and a 2 quantifies the ampli¬ 
tude of the power spectrum relative to the original field. 
Following the same line of reasoning given in the previous 
subsection, we can generalise this result to the anisotropic 
case: 


Pc{k,e) a2P{k,e) . (13) 

As before, the clipped two-point statistics of a second or¬ 
der field maintain the same shape as the undipped case, 
at least for weak transformations. Note that when clip¬ 
ping at equivalent thresholds, uq = the higher order 
field is subject to a significantly stronger suppression of 
its two point statistics than for the Gaussian case, and 
this trend strengthens with yet higher order fields. 


D. Hybrid Fields 

Practical applications of clipping will inevitably in¬ 
volve the superposition of a Gaussian field with other 
components that contaminate the desired signal, partic¬ 
ularly where the amplitude of the field is large. In this 
scenario, clipping can assist in extracting the power spec¬ 
trum associated with the original Gaussian field. 

Consider a hybrid field h{x) that is a linear combina¬ 
tion of a Gaussian field S{x), a higher order field S^{x), 
and a nuisance field 5^ ^ which characterises some un¬ 
known departure from the model: 

h{x) = ai6{x) + a 2 {S‘^{x) — cr^) + S^(x) . (14) 

Upon clipping at a given threshold (5o, and provided the 
nuisance field is constant where h(x) < Sq, we can sub¬ 
tract the mean to remove any residual contribution from 
( 5 ^ 


hc(x) = aiSc + a2(Sc - cr^), 


(15) 


where the component fields Sc and S^ are now also clipped 
fields. This result is important as it shows that we can 
cleanly remove any trace of our nuisance field S ^. In most 
practical applications the recovery will be imperfect, as 
the nuisance field is likely to vary outside of the clipped 
region. However this is a much better state of affairs than 
the conventional approach - be it perturbation theory or 
a model of galaxy bias - where we assume that the extra 
terms missing from our model (as specified by vanish 
everywhere. With clipping, we can now make the much 
more reasonable approximation that S^ only vanishes 
where h(x) is small. 

The two component fields Sc and S^ each experience 
their own distinct thresholds, which may be found by 
solving (15) with the condition hc(x) = Sq. The power 


spectrum of hc(x) is given by 


where Pcl(^) and Pc 22 {k) denote the power spectra asso¬ 
ciated with the clipped S{x) and S‘^{x) fields respectively. 
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The cross spectrum Pci 2 (^) vanishes in the limit of a 
high threshold, but becomes increasingly prominent as 
the threshold is lowered. In order to estimate Pci 2 {k), 
we may decompose it as 

= (n^2) - {r2S^^^) - , (17) 

where the residual fields are defined as Vi = 6^^^ — Sc^. 
Provided the clipping is weak, (5o > cr, the residual fields 
are closely related r 2 — 2 to^ri, where to is the thresh¬ 
old experienced by the Gaussian field. This may be re¬ 
expressed in the form 


Pcl2(k) 




( 18 ) 

Since Pcui^) may be expressed as a linear combination 
of the first and second order power spectra, we can now 
rewrite ([T^ in the form 


Pc(k) aLPhik) + a2P22{k ). 


(19) 


where ul and a 2 are the apparent amplitudes of the orig¬ 
inal linear and second order spectra. This result helps 
explain why the simple model used in Simpson et al 
[9] was particularly successful at reproducing the clipped 
dark matter power spectrum, without explicitly account¬ 
ing for the cross-spectrum Pci 2 {k). 

Adopting this higher order model, rather than relying 
on the linear solution from §IIB holds two advantages. 
First of all weaker clipping thresholds can be used, allow¬ 
ing the power spectrum to maintain a high amplitude. In 
addition, the inclusion of a higher order term potentially 
allows the degeneracy between linear bias and ag to be 
lifted. The disadvantage of this approach is the difficulty 
in estimating a 2 , which could either be calibrated from 
simulations, or simply treated as an additional free pa¬ 
rameter. 


III. GALAXY DENSITY FIELDS 

Galaxy redshift surveys continue to develop an intri¬ 
cate mozaic of the low redshift Universe. By convolving 
the point-like positions of galaxies with a suitable kernel, 
a continuous density field can be generated. These galaxy 
fields are heavily influenced by both nonlinear structure 
and galaxy bias, which have proved highly challenging 
to model. In each case, it is the highest density regions 
which are particularly troublesome, and this motivates 
the application of clipping. The three dimensional na¬ 
ture of the data ensures that clipping can be applied 
very efficiently. Selecting a threshold that affects only 
^ 1 % of the field’s volume typically leads to a reduction 
in large scale power by a factor of two. Maps which are 
two dimensional projections, such as those derived from a 
photometric survey, could also be subject to clipping but 
a greater proportion of the area would need to be clipped 
in order to achieve the same degree of suppression. 


In this section we explore the consequences of clipping 
a galaxy density field in redshift space. We shall work in 
the distant-observer approximation such that all line-of- 
sight displacements may be considered parallel. 


A. Galaxy Bias 

If we model the fractional overdensity of galaxies Sg = 
pglpg — 1 diS an arbitrary function of the local dark matter 
density [23] . 




k=0 


k\ ’ 


( 20 ) 


then applying the clipping transformation suppresses 
higher order terms in the same way that higher or¬ 
der terms in perturbation theory are suppressed. The 
simplest extension to the linear bias model would be the 
introduction of 62 , which for the case of Gaussian dark 
matter fluctuations leaves us with a hybrid field as de¬ 
fined in (14). With a sufficiently low clipping thresh¬ 


old, the linear bias parameter bi dominates such that 
the clipped galaxy power spectrum is highly insensitive 
to the initial value of 62 . This linearisation process was 
demonstrated explicitly in Figure 4 of Simpson et al. [9]. 
Even in the context of more complex models of galaxy 
bias, such as those induced by tidal fields m , we ex¬ 
pect a similar behaviour. Non-linear contributions to the 
galaxy bias still predominantly arise in regions where 5 
is large, and these are the regions suppressed by the clip¬ 
ping transformation. 

There is however a fundamental limit on how much 
we can shield ourselves from the influence of the high¬ 
est density regions, and this stems from the estimation 
of the mean number density. In defining the fractional 
overdensity, 5g = pg/pg — the quantity pg necessarily 
incorporates the abundance of galaxies across the whole 
volume, prior to clipping. Unlike dark matter perturba¬ 
tions where the total particle number is conserved, no 
such restriction applies to galaxy bias. For example if 
baryonic effects reduce the abundance of galaxies in clus¬ 
ters such that the total galaxy count across the survey 
volume is lowered by a small fraction then the inferred 
amplitude of fluctuations across the rest of the volume 
are overestimated by 


K — ^0(1 Py ) 


( 21 ) 


where is the true fractional perturbation. 


B. Redshift Space Distortions: Linear Model 

Whilst the true spatial distribution of galaxies is ex¬ 
pected to be statistically isotropic, their redshift-inferred 
distances receive an additional displacement due to their 
peculiar velocities, generating a statistically anisotropic 
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configuration. This permits a measurement of /ag, 
where the logarithmic linear growth rate is given by 
/ = din J/dlna, and erg defines the amplitude of linear 
perturbations. 

An additional source of anisotropic clustering arises 
from inaccuracies in the assumed geometry of the Uni¬ 
verse, which is required when converting the observed 
values of angles and redshifts into a Euclidean frame¬ 
work. This can potentially generate false measurements 
of the growth rate [25] . In this work we shall consider the 
background expansion to be fixed to a fiat ACDM model 
with Virn = 0-27 unless specified otherwise. 

Clipping in redshift space carries additional compli¬ 
cations. The small scale velocity dispersion associated 
with the ‘Fingers of God’ effect will tend to move galax¬ 
ies out of the high density peaks, and potentially into a 
surrounding region that lies below the clipping thresh¬ 
old. We should therefore expect the removal of non¬ 
linear effects to be less efficient in redshift space. The 
velocity dispersion also causes the power spectrum to 
steepen at larger wavenumbers along the line of sight. 
This strong spectral slope enhances the relative ampli¬ 
tude of the higher order terms in (|^ , so these should not 
be neglected. 

On all but the largest scales, the real space cosmologi¬ 
cal density field at low redshifts is not well described by 
a Gaussian field. It may instead be considered as a su¬ 
perposition of a Gaussian component and an extra 
field 5^ representing the conglomeration of all nonlinear 
corrections. 


sir. +5^. ( 22 ) 

The 5^ field is largest where the linear approximation 
is most strongly violated - both from the truncation of 
higher order terms in perturbation theory and more fun¬ 
damentally from the assumption of a single-valued and 
curl-free velocity field. The matter density field is traced 
by the galaxy number density field, which again may be 
decomposed into a Gaussian component, with a linear 
bias factor 6, and a residual term such that 

= (23) 

Now moving to redshift space, the Gaussian compo¬ 
nent is described by the Kaiser model [26] , which relates 
the real space linear density perturbations with those in 
redshift space, which we couple with a Lorentzian model 
of the velocity dispersion: 




b + fii^ 


yi + fc2(j2^2/2 


Sl^^\k), (24) 


where jj. = k\\/k is defined as the relative fraction of the 
wavevector that extends along the line of sight. The pair¬ 
wise velocity dispersion ay effectively smooths the field 
along the line of sight. Theoretically this quantity is 
given by 


where Pee{k) is the velocity power spectrum [27]. In 
practice we shall treat ay as a free parameter, due to the 
uncertain behaviour of nonlinear motions. 

The clipped galaxy field 5lg can be represented as the 
sum of a clipped Gaussian field and a residual term 5^. 
Therefore the resulting power spectrum may be expressed 
as the sum of three terms, the two autocorrelations and 
the cross-spectrum 


P'lik) = P^^ik) + 2P*i,(fc) + P^{k). (26) 

Given that the nonlinear terms encapsulated by 6^ dom¬ 
inate the clustering statistics at larger values of 6, it ex¬ 
periences a much greater loss of power than the linear 
component. Therefore we should expect that after clip¬ 
ping the first term remains the dominant contribution 
to the total power over a broader range of scales. Our 
simplest model for the clipped galaxy power spectrum is 
therefore encapsulated by 


= C {P^{k,fi);uo} , 


(27) 


where C denotes the transformation defined by Uo 
is the normalised clipping threshold experienced by the 
Gaussian field as given by (|^, and 




2\2 




( 28 ) 


where Pg{k^ /i) is the Gaussian contribution to the galaxy 
power spectrum in redshift space, Pi,{k) is the linear mat¬ 
ter power spectrum in real space, and the anisotropy pa¬ 
rameter /3 = f /b quantifies the level of anisotropy in the 
galaxy power spectrum induced by linear velocity pertur¬ 
bations. 

For the case of a hybrid field such as the one defined 
in (14), then as clipping is applied, and nonlinear con¬ 


taminations are suppressed relative to the linear contri¬ 
butions, we should expect the recovered value of (3 to 
be closer to the theoretical value for a given /Cmax; alter¬ 
natively, we should be able to achieve the same level of 
systematic error in at a higher /cmax- The actual level 
of error and/or smallest scale to probe must be deter¬ 
mined empirically using simulations, as described in the 
following section. 


C. Redshift Space Distortions: Nonlinear Model 

Our base model is defined by the set of four parameters 
6(78, and relies upon the recovery of the linear 

matter power spectrum. However the linear power spec¬ 
trum decays rapidly towards higher wavenumbers, and 
by k r\j 0.5 hMpc ^ is typically an order of magnitude 
lower than the contribution from the one-loop correction 
to the power spectrum. Therefore despite the suppres¬ 
sion of higher order terms, the inclusion of a suppressed 
one-loop term substantially improves the model for the 
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real space power spectra of matter and galaxies [9] . Fur¬ 
ther motivated by the results of we introduce an 

extended model with the additional parameter e, which 
accounts for a higher order contribution to the power 
spectrum 


b‘^(l + 

= °^ l + (fcaS)V2 + e^’iioop(fc)] . (29) 


The parameter ul is the coefficient of the linear power 
which contributes to the clipped power spectrum, as de¬ 
fined in the model of (19). In general the value of ul 


cannot be determined a priori. It may be evaluated em¬ 
pirically by considering the fractional change in ampli¬ 
tude of the large scale clustering of the field, 


frequencies of = 0.8hMpc“^ and 1.6hMpc“^ respec¬ 
tively. The standard deviations are 3.1 for the coarse grid 
and 6.0 for the high resolution field. 

In order to evaluate the true value of (3 we need the 
linear bias parameter. The ratio of the real space galaxy 
power spectrum with that of the Millennium simulation’s 
dark matter field gives h = 1.125 ± 0.01 at = 0, on the 
largest available scales. At 2 : = 0.687 this increases to 6 = 
1.29 ± 0.02, while the growth rate is well approximated 
by f{z) z, 


B. Threshold Selection 


Ul 


lim 

r^oo 


Ur) 
U) ' 


(30) 


IV. SIMULATIONS 

In order to test our theoretical models, we construct 
mock density fields for both dark matter and galaxies. 
These are derived from numerical simulations, and trans¬ 
formed into redshift space using the distant observer ap¬ 
proximation. In this section we summarise our methods 
for generating and modelling the power spectra associ¬ 
ated with clipped cosmological fields, and for the estima¬ 
tion of their covariance matrices. 


A. Number Density Fields in Redshift Space 

For our mock dark matter field we utilise the z = 0 
snapshot from the Horizon Run 2 simulation [28], which 
consists of 6,000^ particles within a periodic box of size 
7,200 h“^Mpc. The amplitude of linear perturbations is 
(78 = 1/1.26, with a matter density ftrn = 0.26. The 
redshift-space density field is defined by considering a 
1920^ grid using a Nearest Grid Point (NGP) scheme, 
where the particles are displaced along one axis in accor¬ 
dance with their peculiar velocity. This leaves ns with 
a grid size of 3.75h“^Mpc corresponding to a Nyquist 
frequency of k-^ = 0.84hMpc“^. 

The mock galaxy catalogues are taken from Guo et al 
[29] . which applies a semi-analytic model to the halo 
merger trees of the Millennium-I simulation [30] . Follow¬ 
ing Simpson et al. m we use the two snapshots at ^ = 0 
and = 0.687 in order to explore different amplitudes of 
matter perturbations and growth rates. With ag = 0.9, 
the 2 : = 0 snapshot possesses a slightly higher amplitude 
of fluctuations than the Horizon simulation. Applying a 
stellar mass cut of log;Lo(^*/^0^) ^ ^ leaves ns with a 
distribution resembling that found in our GAMA sam¬ 
ple [31]. Number density fields are formed on both 128^ 
and 256^ grids across the 500 h“^Mpc box, with Nyquist 


After constructing each number density field, we apply 
the clipping transformation defined by 0. In order to 
select a suitable threshold value ^0 we require an appro¬ 
priate metric for defining the strength of clipping. For 
the case of a Gaussian field, the normalised threshold 
1^0 = (5 o/a/2(7 provides a natural measure for this in re¬ 
lation to the standard deviation of the field. However 
when working with fields that are highly non-Gaussian, 
we do not know a priori what impact a given threshold 
will have. Some degree of iteration is therefore required 
in order to reach the desired reduction in power, charac¬ 
terised by ul- 


When working with the simulations, rather than quot¬ 
ing absolute values of the threshold (5o, we choose to work 
in terms of the fraction of mass (or galaxies) removed. 
This way fields with larger fluctuations naturally adopt 
higher thresholds (5o- 

The disadvantage of a stronger (lower) threshold is 
a larger drop in the amplitude of the power spectrum, 
which in turn reduces the maximum wavenumber avail¬ 
able before the shot noise contributions dominate. Very 
strong thresholds also induce a large contribution from 
the cross-power Pi 2 {k), as given by (18), further re¬ 
ducing the power. The optimal choice of threshold is 
therefore one that adequately suppresses contributions to 
the power from nonlinear structure and nonlinear bias, 
while maintaining a relatively high amplitude of the lin¬ 
ear power spectrum. As shown in Figures 6, 7, and 8 of 
Simpson et al. m, the decay of the higher order term 
rapidly outpaces the contribution from the linear power 
spectrum, such that a factor of two reduction in linear 
power is sufficient to eliminate approximately 80% of the 
nonlinear power. 


For each field we explore a range of threshold values, 
since it is important to verify that different thresholds 
generate consistent parameter constraints. Performing a 
likelihood analysis that combines the power spectra de¬ 
rived from different thresholds may further improve pa¬ 
rameter constraints. However an estimation of the covari¬ 
ance between the different spectra is beyond the scope of 
this work, here we shall only consider the analysis of each 
threshold separately. 
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C. Methods 


The theoretical model for the clipped power spectrum 
is evaluated with the following procedure 


Apply the redshift space model defined in (28) to 
determine /i) from the linear power spectrum 

evaluated by CAMB [32] combined with the chosen 
parameter values of 6, /, and Gy. 


• Interpolate the power spectrum onto a 3D grid 
matching the specifications of the NGP lattice gen¬ 
erated from the simulations. 


• Transform to the 3D correlation function ^(cr, tt) 

• Evaluate the clipped 3D correlation function 
^c(cr, tt) using We find it sufficient to truncate 
the series expansion at rimax = 4 in order to estab¬ 
lish sub-percent accuracy. This step is discussed 
further in Appendix [A| 

• Apply the inverse Fourier transform to determine 
the expected 3D power spectrum as a function 
of wavenumber magnitude k and orientation 

The power spectrum is assigned to lin¬ 
early spaced bins, with widths of Ak = 0.025, 
A/i = 0.2. We maintain the same binning 

scheme throughout this work. Note that the mean 
wavenumber contributing to the bin is in general 
at higher values of k than the bin centre, due to 
the abundance of modes within a spherical shell dk 
scaling as k‘^. 


A four-dimensional likelihood grid is constructed, 
jC{Pc{k, /i) jttL, /d, bag)^ where the parameter set relates 
to the model given by (27) and (28). 

With the clipped power spectrum alone as the only 
source of information, we would have little knowledge of 
how strong the applied clipping has been, and therefore 
ul is poorly constrained. However the drop in power 
relative to the original (unclipped) power spectrum pro¬ 
vides us with some extra information that can assist 
in constraining the range of ul values. We therefore 
make use of the fractional loss of power experienced at 
k < O.lhMpc”^ and fi < 0.6, to provide some external 
information for the value of ul- At higher values of k 
and fi the larger contributions from nonlinear structure 
mean that they experience a greater drop in power. 

We assign flat priors of 0 < ul < 1, 0 < /3 < 2, 0 < 
hag < 2, and for the velocity dispersion we use a broad 
Gaussian prior = 300 ± 300 kms“^, as defined in 

(28). Our results are largely insensitive to the particular 
choice of priors. 


D. Covariance Estimation 

The large volume of the dark matter field in the Hori¬ 
zon Run 2 simulation, (7,200 h“^Mpc)^, permits a direct 


estimation of the 70 x 70 covariance matrix associated 
with the P(/c,/i) bins by considering the covariance of 
power spectra evaluated from 350 subvolumes, each of 
size (960 h“^Mpc)^. To evaluate the covariance matrix 
we employ the Ledoit-Wolf shrinkage estimator [33], fol¬ 
lowing the prescription 

E = (i^F + (l-(i^)S, (31) 

where F is our estimated covariance matrix, and is 
the shrinkage constant. S is the sample covariance ma¬ 
trix, defined as the ensemble average over N = 350 sub¬ 
volumes 


Si, = {APiAPj ), (32) 

where APi represents the deviation in power in the ith 
P{k, fi) bin from the sample mean. The shrinkage target 
F is defined in terms of the sample covariance 


Fa = < 


Sij , 


{i=J) 


{i ^ j) 


(33) 


r = 


AT-l N ^ 


(34) 


where N = 350. The shrinkage constant is estimated 
using publicly available code m- We emphasise that 
this approach yields a considerably improved estimate 
of the covariance matrix than simply using the sample 
covariance S alone. 

For the galaxy sample associated with the Millennium 
simulation, the volume is considerably smaller and there¬ 
fore a direct estimate of the covariance matrix would be 
prohibitively noisy. Instead we only estimate the diago¬ 
nal terms of the covariance matrix, again using the vari¬ 
ance of sub volumes. The estimation of off-diagonal terms 
is determined by using the dark matter covariance as a 
template, such that 

Cij = coTI{^,j)^yCiiCjj . (35) 

where corr(i,jf) is the correlation matrix from the dark 
matter power spectrum. 


E. Shot Noise Estimation 

For the conventional unclipped power spectrum, the 
discrete nature of sources leads to an additional contri¬ 
bution of power, which if we assume to be Poissonian in 
nature is given by Pg = n~^. As clipping smooths the 
field above the threshold value, the shot noise contribu¬ 
tion is reduced. Approximating the noise field as Gaus¬ 
sian allows us to utilise ^ to estimate Pg — 
where fv is the fraction of the volume of the field lying 











below the clipping threshold. In this work the typical 
volume fraction is of the order ^1%, and therefore the 
correction to the shot noise is negligible. 

The power spectra of clipped fields are highly robust 
to changes in the number density of sampled points [9]. 
The only noticeable consequence appears to be that fields 
with higher shot noise possess noisier power spectra and 
can therefore not utilise as wide a range of wavenumbers. 


V. RESULTS FROM SIMULATIONS 

First we present the power spectra from the dark mat¬ 
ter and galaxy fields, at different clipping thresholds, be¬ 
fore reviewing the results of the likelihood analysis. 


A. Clipped Power Spectra 


The panels in the left hand column of Figurecompare 
the dark matter power spectrum from the Horizon Run 
2 simulation with the linear theory prediction in each of 
five angular bins. Within each panel the uppermost set 
of points represent the power from the original undipped 
field, while the central set is generated after applying a 
logarithmic transformation ln(l + S). The lowest set of 
point corresponds to a field subject to a clipping trans¬ 
formations, Q, with a threshold value (5o chosen such 
that 20% of the mass is removed. Each dashed line cor¬ 
responds to the model based on the linear power spec¬ 
trum given by (28) where the amplitude is rescaled to fit 


the transformed spectra. The form of the real space lin¬ 
ear power spectrum and the anisotropy parameter /3 are 
assigned values according to the simulation parameters. 
The value of P = f /b is derived using the linear growth 
rate / and the linear bias 6=1 since we are 

working directly with the dark matter. No error bars are 
displayed in these panels because the statistical error is 
considerably smaller than the marker size. The bins in 
k and fi were chosen to match the power spectra derived 
from the GAMA survey. 

For such an evolved field the linear theory prescription 


given by (28) typically holds only on very large scales. 


Even at k O.lhMpc” the model overestimates the 
power in the highest /i bin by almost 10%, consistent with 
the findings of Jennings et al. [35]. However once either 


transformation is applied, the linear formalism of (28) 


provides a significantly improved description. Agreement 
with the clipped spectrum is better than 5% within the 
range {k < 0.5hMpc“^;/i < 0.8). This improvement 
in the modelling occurs due to the strong suppression 
of higher order terms in perturbation theory |9|. The 
leading cause of tension with the model appears to be 
within the highest /i bin, which is perhaps unsurprising 
since these modes receive contributions from very small 
physical scales, due to the velocity dispersion of galaxies. 
The central set of points illustrate the response of the 
power spectrum to another local transformation, the log¬ 


arithm of the number density, ln(l + 5). Neyrinck et al 
|8| demonstrated that this can help linearise the power 
spectrum of the real space dark matter field. We find that 
considerable linearisation also occurs when applying the 
log transform to the dark matter field in redshift space. 
The shape of the linear theory power spectrum defined by 
(28) is recovered to better than 10% for k < 0.4hMpc“^. 

The central column of Eigurel^is in the same format 
as the left column, but illustrates the galaxy power spec¬ 
trum from the Millennium Simulation at ^ = 0 before 
and after applying transformations to the number den¬ 
sity field. As before, the uppermost set of points in each 
panel represent the original unclipped field. The lowest 
set of points are generated by clipping 20% of the galax¬ 
ies, which brings the data points closer to the shape of 
the linear model (28), as given by the dashed line. To 
obtain the value of /3 for the model requires a combina¬ 
tion of the linear bias, which is estimated from the ampli¬ 
tudes of the largest Eourier modes in the simulation to be 
6 = 1.125, and the growth rate f{z = 0) = 0.47. As quan¬ 
tified in MV A[ the true value of the linear bias is only 
an estimate, however it remains a subdominant source 
of uncertainty. The error bars in the central and right 
hand columns are significantly larger than those in the 
left hand column, reflecting the considerably smaller box 
of the Millennium-I Simulation (500 Mpc/h) compared 
to that of the Horizon Run (7,200 Mpc/h). 

Deviations between the clipped spectrum and the lin¬ 
ear model remain lower than 10% for all data points at 
k < 0.3hMpc“^. Unlike the case of dark matter, it is 
the lowest /r bin that causes the greatest tension with 
the model. This may be due to the local motions of 
galaxies causing them to be displaced from their high 
density regions, which would make the clipping process 
less efficient, leaving behind a considerable proportion of 
the nonlinear contributions to the power spectrum. The 
middle set of points in Eigure correspond to the loga¬ 
rithmic transformation, but now only a modest degree of 
linearisation is observed. This reduced effectiveness can 
be attributed to the sampling noise from the galaxy fields 
with n = O.OlMpc”^. Due to this sensitivity to the level 
of shot noise, we shall focus on the clipping transform for 
the remainder of this work, which by contrast is largely 
insensitive to shot noise. 

The right hand column of Eigure[^ explores the power 
spectra for a different sample of galaxies at a different 
redshift, 2 : = 0.687. Here both the linear bias 6 = 1.29 
and growth rate /(z = 0.687) = 0.77 have changed from 
those of the central column, yet the outcome is similar, 
in that the application of clipping significantly improves 
the performance of the Kaiser model defined in (28). As 
with the low redshift galaxy sample, the departure from 
linearity is less than 10% for k < 0.3hMpc“^. 

Eigurej^ displays three subsections of the 70 x 70 cor¬ 
relation matrix (14 bins in /c, 5 bins in /i) associated with 
the unclipped dark matter field. The coupling of modes 
becomes particularly apparent towards higher wavenum¬ 
bers, k > 0.25hMpc~^. It is interesting to note that 
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FIG. 1. The redshift space power spectra P{k,/j) from the simulations described in Section [Tv| divided into five angular bins. 
The three columns are dark matter from the Horizon simulation (left); galaxies at z = 0 from the Millennium simulation 
(middle); and galaxies at z = 0.687 from the Millennium simulation (right). Within each panel the original power spectrum 
is shown as black set of data points while the red and blue sets correspond to the power spectra generated after the number 
density field has been subject to logarithmic and clipping transformations respectively. The clipping threshold is chosen such 
that it lies below 20% of the mass. The dashed lines represent the linear model specified by (28). For the case of the transformed 
spectra, the amplitude of the model is rescaled to fit the data. 
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FIG. 2. Three subsections of the 70 parameter correlation matrix from the power spectrum P(/c, ji) of the undipped dark matter 
field in the Horizon simulation. The left hand panel corresponds to the lowest bin in wavevector orientation, 0 < /x < 0.2, 
while the centre panels relates to the highest bin, 0.8 < /x < 1. Each section consists of ten bins in wavenumber k. The right 
hand panel is subsection of the correlation matrix associated with the cross-correlation between those two /x bins. The growing 
off-diagonal correlations towards smaller scales, particularly apparent in the lowest /x bin, is an indicator of the nonlinear 
structure that generates mode coupling. The correlation matrix was estimated by following the prescription described in §IV D[ 


these off-diagonal terms fall by approximately 20 — 30% 
after clipping has been applied. This decorrelation of 
neighbouring bins was previously observed in the real 
space power spectrum [9]. 

The shrinkage constant also reduces considerably after 
the application of clipping. From the undipped field we 
find a shrinkage constant of 5'^ = 0.25. Thresholds se¬ 
lected to remove 5% and 10% of the dark matter yields 
shrinkage constants of 0.18 and 0.14 respectively, reflect¬ 
ing the increasingly Gaussian nature of these fields. 

In Figure we explore the efficacy of fitting the power 
spectra of clipped fields with the non-linear model, spec¬ 
ified by ( [^ , which is based on the model of the real 
space power spectrum presented in Simpson et al. [9]. 
The simulations in question are the same as those used 
in Figure The data points used to fit the model pa¬ 
rameters are k < 0AhMpc~^ and fi < 0.8. The upper 
set of data points corresponds to the original field, with 
the dashed line now making use of the one-loop power 
spectrum. The middle and lower sets of data points re¬ 
late to clipping 10% and 20% of the mass respectively. 
Now that an extra contribution from the one-loop power 
spectrum is included, the data points in the lower /x bins 
are much better accounted for, compared with the linear 
model in Figure However the highest /x bin appears 
significantly underestimated. 

We conducted further investigations by evaluating the 
power spectra of dark matter haloes in the Millennium 
simulation, before and after clipping, at the same two 
snapshots as the aforementioned galaxy catalogues. A 
very similar trend is found, whereby the clipped spectra 
are readily described by the model for k < 0.4/iMpc“^. 
Since these power spectra appear very similar to those 
displayed in Figure they are not shown here. 


B. Growth of Structure 

First we attempt to recover the correct cosmological 
parameters from the clipped power spectra, using only 
the linear model given by ( [^ : 

(36) 

Our basic set of parameters is {/5, cr^, ap}. The nor¬ 
malised growth rate fag may be derived from these via 
the following relation: 

fas = pbas ■ (37) 

Figure]^ shows the constraints on fas when applying the 
linear model to the three different fields from the simula¬ 
tions, as a function of the maximum wavenumber /Cmax- 
In each case the true value of fag is illustrated by a hori¬ 
zontal dashed line. The points have small horizontal off¬ 
sets for clarity, and appear in order of increasing clipping 
strength. The squares, triangles and inverted triangles 
correspond to clipping thresholds below 5, 10 and 15% of 
the field respectively. The shaded regions represent their 
68% confidence limits. For reference, the maximum like¬ 
lihood points from the original unclipped field are shown 
as black circles. 

Results from the dark matter field are shown in the 
left hand panel. Removing only 5% of the mass is found 
to be sufficient to correct for much of the nonlinear be¬ 
haviour on scales k < 0.3/iMpc~^. Similar behaviour 
is seen in the central panel of Figure which uses the 
2 ) = 0 galaxy sample, with a true value of erg = 0.9. In 
the right hand panel we find that the galaxy sample at 
^ = 0.687 also significantly improves the recovery of the 
underlying cosmology when using the linear model on 
scales k < 0.3/iMpc“^. The shaded regions are signifi¬ 
cantly broader in the central and right hand panels, re¬ 
flecting the considerably smaller box of the Millennium-I 
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FIG. 3. A similar format to Figure except here the dashed lines represent the best-fit non-linear model, as given by (29). 
The central and lower sets of data points correspond to different clipping thresholds, removing 10% and 20% of the mass 
respectively. 
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FIG. 4. Marginalised constraints on the normalised growth rate fag from the simulations, when using only the linear power 
spectrum \21\ to model the data. The three shaded regions arise from the power spectra after clipping 5% (squares), 10% 
(triangles), and 15% (inverted triangles) of the mass from each held. The horizontal dashed lines indicates the true values of 
fag within each snapshot. 



FIG. 5. The same format as Figure except we now ht the data using the 
prior — 1 < e < 1. 


nonlinear model dehned by (29), with a uniform 


Simulation (500^Mpc) compared to that of the Hori- C. Nonlinear Model 

zon Run (7, 200/i“^Mpc). 


In the context of the more general model ( [2^ , which 
invokes additional contribution controlled by the e pa¬ 
rameter, we know that for very weak clipping e will be 
positive and when the clipping is very strong e becomes 
negative as the contributions from cross spectra such as 
Pi 2 {k) dominate. Therefore there is inevitably a thresh¬ 
old at which e vanishes and the linear model offers a good 
fit to the data. This is represented by the inverted trian¬ 
gles in Figure]^ 


We repeat our analysis using the extended model de¬ 
fined by (29), which introduces an extra parameter e to 
control the amplitude of the clipped one-loop power spec¬ 
trum, which is otherwise fully specified in terms of the 
linear power spectrum. Figure shows the constraints 
derived on fag from the three simulated fields clipped 
with the same thresholds as Figure The extra free¬ 
dom in the real space matter power spectrum leads to 
significantly improved measurements at weaker clipping 
thresholds. In the case of dark matter we find that the 
maximum likelihood is within 10% of the true value for 
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FIG. 6. The same format as Figure [^except here we only make use of intermediate wavevector angles (0.2 < fi < 0.8). 


each /cmax < 0.4/zMpc“^ and for each threshold. Simi¬ 
larly in the galaxy field both clipping thresholds return 
more consistent constraints, and with only a modest loss 
of precision compared to the simpler linear model. How¬ 
ever we find that the extreme values of /i are responsible 
for the bulk of the tension between the model and the 
data. It may be the case that more complex models such 
as those outlined by Taruya et al [36] may provide a bet¬ 
ter description of the anisotropies in the clipped power 
spectrum. An exploration of these models in the context 
of clipped fields is beyond the scope of this work. Re¬ 
stricting ourselves to 0.2 < /i < 0.8 improves agreement 
between the model and data to better than 2% across 
all wavenumbers k < 0.5/iMpc~^. This is reflected in 
Figure where the tendency to overpredict fag at the 
highest values of /cmax is resolved. 


VI. DATA 

The galaxy redshift surveys that have the greatest po¬ 
tential to benefit from clipping are those with a high num¬ 
ber density of galaxies. A densely sampled field ensures 
that shot noise is low out to high wave numbers, even 
after the drop in the amplitude of the power spectrum 
due to clipping. The GAMA survey provides an excellent 
basis for the first application of clipping to a real galaxy 
field. In this section we present details of the survey, and 
how the power spectra were generated. 


A. The GAMA Survey 

The Galaxy and Mass Assembly (GAMA) project [37b 
l4Q] is a multi-wavelength photometric and spectroscopic 
survey. The redshift survey, which has been carried out 
with the Anglo-Australian Telescope (AAT), has pro¬ 


vided a dense, highly-complete sampling of large-scale 
structure up to redshift ^ ^ 0.5. The primary target 
selection is r < 19.8 (where r is an extinction-corrected 
SDSS Petrosian magnitude), using TilingCatv41. 

Following Blake et al m, we analyzed a highly- 
complete subsample of the survey dataset known as the 
GAMA II equatorial fields. This subsample covers three 
12 X 5 deg regions centred at 09h, 12h and 14h30m, which 
we refer to as G09, G12 and G15, respectively. Galaxy 
redshifts were obtained from the AAT spectra using a 
fully automatic cross-correlation code that can robustly 
measure absorption and emission line redshifts [42]. We 
restricted the redshift catalogue to galaxies with “good” 
redshifts (NQ > 3). In order to obtain high-resolution 
measurements of the density field we restricted our anal¬ 
ysis to the redshift range 0.002 < z < 0.25, where 
the galaxy number density exceeds 10“^ Mpc“^. In 
the (009, G12, G15) regions we utilized (32076, 37382, 
36538) galaxies in our analysis. The comoving volume of 
each region is approximately 6.4 x 10^(/i“^Mpc)^ [41]. 

The survey selection function at each point, used in the 
calculation of the galaxy overdensity, was determined by 
combining the angular completeness map of the survey 
(which has a mean value of 97% across the three regions) 
with an empirical fit to the galaxy redshift distribution, 
performed after stacking together the data in the three 
regions to reduce fluctuations due to cosmic variance. 
Full details of the method are described in Section 3.2 of 
Blake et al [41]. 

B. Estimating the Clipped Power Spectrum 

The clipped power spectra for each GAMA region were 
determined for a given overdensity threshold 6c\ip as fol¬ 
lows: 

1. The galaxy distribution was binned on a common 
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3D grid to the selection function, with a resolu¬ 
tion of 2h~^ Mpc. We denote the gridded distri¬ 
butions from the data and random samples as D{x) 
and R{x)^ respectively. The random catalogues are 
sampled from the selection function constructed for 
the GAMA survey data, which combines the an¬ 
gular completeness in each survey region with an 
empirical smooth redshift distribution fit to a com¬ 
bination of the three regions. 

2. The distributions were smoothed using a Gaussian 

kernel G{x) = ^ yVe take A = 2h~^ Mpc 

for our analysis. We denote the smoothed fields by 
sm{D) and sm{R). 

3. The overdensity field for each region was estimated 
as S{x) = sm{D)/sm{R) — 1, where the normaliza¬ 
tion of R was fixed such that {6) =0. 

4. The mean overdensity of each region in the redshift 
range z < 0.25, relative to the average of all three 
regions, was estimated using the measured number 
of galaxies as (5reg = (—0.0922,0.0580,0.0341) for 
(GOO, G12, G15). The effective clipping threshold 
applied to the locally defined fluctuations within 
each region was then adjusted to (5ciip,eff to ac¬ 
commodate these mean density fluctuations, where 

1 T ^clip,eff — (t T ^clip)/(f T ^reg)- 

5. For any grid cell with S > ^ciip,eff: the unsmoothed 
gridded data value was lowered to D = R{1 -\- 

^clip,eff ) • 

6. The power spectrum of the clipped gridded data 
held D was measured using Fast Fourier Trans¬ 
form (FFT) techniques following Section 3.3 of 
Blake et al. m- The optimal-weighting estimation 
scheme of Feldman, Kaiser & Peacock (1994) was 
applied, assuming a characteristic power spectrum 
amplitude P{k) = 5000 Mpc^. We binned the 
power spectrum by k and /i, where /i is the cosine 
of the angle of the wavevector k with respect to the 
line-of-sight, using bin widths Ak = 0.05 h Mpc“^ 
and A/i = 0.2. The integral constraint correction 
to the power spectrum was included in the esti¬ 
mation process (using the Fourier transform of the 
window function). 

7. The amplitude of the measured power spectra was 
corrected for the mis-estimate of the mean density 
of the region, through multiplication by a factor 

(1 + ^reg)^- 

8. The convolution matrix, which is used to project a 
model power spectrum P(/c, /i) to form a compar¬ 
ison with the data given the survey window func¬ 
tion, was determined using the method outlined in 
Section 3.3 of Blake et al [41], in which the full FFT 
convolution is applied to a series of unit model vec¬ 
tors, and an equivalent matrix is constructed row- 
by-row. 


9. The covariance matrix of the power spectrum mea¬ 
surement in (/c,/i) bins was estimated by evaluat¬ 
ing the sums described in Section 3.4 of Blake et al. 
m- Initially the measured power spectrum in each 
bin was used to specify the cosmic variance compo¬ 
nent. This produces an error estimate that is corre¬ 
lated with the data. To resolve this we modified the 
computation using an iterative procedure in which 
the best-fitting (convolved) theoretical model was 
determined and the covariance was re-estimated us¬ 
ing that model. Two iterations were used to ensure 
convergence. 

We repeated the above analysis for three different clip¬ 
ping thresholds (5o = {8:5,4}. These values were se¬ 
lected on the basis of generating a suppression of linear 
power between 30 and 60%. This provides an appropri¬ 
ate balance between the elimination of nonlinear struc¬ 
ture and maintaining a high degree of signal to noise. 
The three thresholds affect approximately 0.7%, 2.1%, 
and 3.2% of the held in terms of volume, and approx¬ 
imately 6.4%, 15%, and 20% in terms of galaxies. The 
latter quantity is defined as the fractional reduction in 
the value of D{x) due to clipping. 

Unlike non-local transforms such as those used for re¬ 
constructing the baryon acoustic oscillations, clipping 
commutes with the window function. This facilitates our 
interpretation of clipped power spectra, since the clipping 
transformation associated with the full (non-windowed) 
universe can be evaluated first, before compensating for 
the impact associated with the window function of the 
survey. 

VII. RESULTS FROM GAMA 

In this section we perform a likelihood analysis to es¬ 
timate the normalised growth rate fag at 2 ; = 0.18. The 
effective redshift of power spectrum measurements in the 
GAMA regions was determined by Blake et al. [41]. The 
methodology from Section |IV| is applied to each of the 
clipped galaxy power spectra from each of the three fields 
of the GAMA survey. 

Figure illustrates the effect of the clipping transfor¬ 
mation. The solid bars represent the probability density 
function of the galaxy density held within the G09 region, 
defined in terms of 40 bins which are spaced equally in 
log(l + J). The hollow bars show the resulting distribu¬ 
tion function if we then apply a clipping transformation 
with a threshold of (io = 4. The distribution function 
remains unaltered below the threshold value while all 
contributions from greater overdensities are compressed 
into the bin associated with the threshold value. 


A. Clipped Power Spectra 

The panels in Figure [^illustrate the anisotropic power 
spectra derived from the three fields (G09, G12, and 
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FIG. 7. The solid bars represent the probability density func¬ 
tion for the galaxy density field within the G09 GAMA region. 
The hollow bars demonstrate the effect that applying a clip¬ 
ping threshold of 5o = 4 would have on this field. A vertical 
dotted line demarcates the underdense and overdense regimes. 
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TABLE I. The maximum likelihood fag values and their 
associated 68% confidence intervals, when fitting the linear 
model to the GAMA data. Each column represents a differ¬ 
ent clipping threshold while each row represents a different 
kmax condition. Table elements associated with a value of 
Xred <1-2 are highlighted in bold. There is a strong covari¬ 
ance between nearby table elements, so we do not attempt to 
combine constraints from different clipping strengths in this 
work. 


G15). As with the simulations, the power spectrum is di¬ 
vided into five equal bins in /i, spanning 0 < /i < 1, while 
the wavenumber bin width is taken to be Ak = 0.05. 
Within each field the three sets of points correspond to 
the power in the field before (black) and after the ap¬ 
plication of clipping thresholds Jq = 8 (red) and Sq = 4 
(blue). At each clipping strength, the dashed line reflects 
the linear model, with the maximum likelihood values of 
/3, cr^, and bag. Estimates of the ap parameter are shown 
in Figure These are determined by the fractional drop 
in P{k, fi) for k < 0.1 hMpc~^ and /i < 0.6, after clipping 
is applied. 


B. Linear Model 


^max 

Undipped 

00 

II 

II 

II 

0.1 

2.77 

1.50 

1.31 

1.26 

0.2 

3.68 

1.00 

0.88 

0.73 

0.3 

2.81 

1.26 

1.05 

1.08 

0.4 

2.66 

1.35 

1.01 

1.00 

0.5 

3.06 

1.87 

1.31 

1.21 


TABLE 11. The reduced values for a variety of clipping 
thresholds and maximum wavenumbers, when fitting the lin¬ 
ear model to the GAMA data. 


0.18). Here we shall present results which combine the 
likelihoods of the three regions. Individual results from 
the three separate regions can be found in Appendix [C] 


First we employ the linear model defined by and 


(27), while fixing the shape of the linear power spectrum 
to the fiducial model. The left hand panel of Figure [To] 
show the maximum likelihood values and 68% error bars 
associated with fag, under a range of different clipping 
thresholds and k^ax values. As before, the squares, trian¬ 
gles and inverted triangles correspond to clipping thresh¬ 
olds of (5o = 8, 5 and 4 respectively. The shaded regions 
represent their 68% confidence limits. For reference, the 
maximum likelihood points from the original unclipped 
field are shown as black circles. Their confidence lim¬ 
its are suppressed for clarity, as they do not provide an 
acceptable fit to the data. 


With the unclipped data the constraint on fag is 
highly sensitive to variations in k^ax, which is consistent 
with the behaviour found in the simulations. Since the 
model inevitably underestimates the amount of real space 
power towards larger k, this leads to an under-estimation 
of ay which in turn biases the estimate of P to be low. 
Applying a high clipping threshold ((^o = 8) show a mod¬ 
est improvement in terms of consistency and goodness of 
fit. Stronger thresholds of Jo = 5 and Sq = 4 provide a 
much improved agreement with the model, and more con¬ 
sistent results towards higher wavenumbers. The maxi¬ 
mum likelihood values are displayed in Table [Tj and the 
reduced values are can be found in Table [ hJ The power 
spectra associated with the clipped fields are found to ad¬ 
here to the linear theory prediction more closely than the 
original field. From the simulations we expect a signif¬ 
icant systematic error to arise at k^^x ^ 0.3/iMpc~^, 
so for the linear model we use k < 0.2hMpc~^ to find 
fag{z = 0.18) = 0.26 ± 0.10. Constraints from the 
clipped fields show a more consistent result across the 
range of wavenumbers than the original field, and also 
have a much improved goodness of fit. 


Following the procedure outlined in Section |VB| we 
use the clipped power spectra from the three GAMA re¬ 
gions to measure the normalised growth rate fag{z = 


Extracting robust constraints from higher wavenum¬ 
bers requires a higher order model, since at these scales 
the amplitude of the linear power spectrum falls far below 
the non-linear contributions. 
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FIG. 8. The galaxy power spectra Pg(k^ii) measured within the three GAMA regions, each divided into five angular bins. 
Within each panel the original power spectrum is shown as the upper set of data points while the middle and lower sets 
correspond to the power spectra generated after the number density field has been clipped with thresholds of (5o = 8 and (5o = 4 
respectively. These thresholds generate a drop in large scale power of approximately 30% and 60%. The upper dashed line 
is generated from the conventional model given by (28). The lower dashed line is the fiducial model for the clipped linear 
power (27). At each clipping strength [3 and ay are set to their maximum likelihood values. The theoretical lines differ slightly 
between each region since each has been convolved with the window function of their respective region. For visual clarity a 
small horizontal offset is applied to the black and blue data points. 
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1 


0.8 


0.6 

L 

0.4 


0.2 


0 



009 012 015 


^max 

Undipped 

5o 


: 8 

<5o 


5 

5o 


4 

0.1 

0.41 

±0.07 

0.34 

± 

0.19 

0.35 

± 

0.19 

0.32 

± 

0.19 

0.2 

0.39 

±0.08 

0.23 

± 

0.11 

0.27 

± 

0.11 

0.25 

± 

0.11 

0.3 

0.27 

±0.07 

0.27 

± 

0.11 

0.30 

± 

0.11 

0.29 

± 

0.10 

0.4 

0.28 

±0.05 

0.20 

± 

0.08 

0.22 

± 

0.08 

0.21 

± 

0.07 

0.5 

0.34 

±0.04 

0.32 

± 

0.08 

0.32 

± 

0.10 

0.30 

± 

0.09 


TABLE III. The same format as 
nonlinear model defined by (29). 


Table but now using the 


^max 

Undipped 

00 

II 

II 

II 

0.1 

2.77 

1.50 

1.31 

1.26 

0.2 

3.12 

0.96 

0.88 

0.76 

0.3 

2.04 

1.36 

1.12 

1.12 

0.4 

2.34 

1.39 

1.02 

1.02 

0.5 

2.55 

1.86 

1.42 

1.27 


FIG. 9. Estimates of the fractional drop in linear power in 
each of the three GAMA regions as a result of applying clip¬ 
ping thresholds of (5o = 8 (upper points) and — A (lower 
points). 


C. Non-linear Model 


The central panel of Figure shows the constraints 
on fag when using the extended model defined by 
(29). Again we find consistent behaviour between dif¬ 


ferent clipping strengths, and across a variety differ¬ 
ent maximum wavenumbers. The extra degree of free¬ 
dom does not appear to significantly weaken the con¬ 
straints. Guided both by the performance of simula¬ 
tions, and the goodness of fit between the data and 
the model, we adopt our benchmark measurement to be 
fcFs{z = 0.18) = 0.29 ± 0.10, using k < 0.3/iMpc“^ and 
So = 4. This measurement is consistent with that derived 
from the linear model, and serves as the central result of 
this work. Our result is consistent with the findings of 
Blake et al. [41], who used the same (unclipped) galaxy 
field to determine f(Ts{z = 0.18) = 0.36 ± 0.09. The full 
set of constraints on fag is presented in Tab le |III[ while 
the reduced values are displayed in Table [iVp^ 

A more conservative approach is to restrict our analy¬ 
sis to intermediate wavevectors, 0.2 < /i < 0.8, and the 
results are shown in the right hand panel of Figure 
While the susceptiblility to systematic errors has been re¬ 
duced, there is also a substantial loss of precision. There¬ 
fore even when making use of the full range of wavenum¬ 
bers, k < 0.5/iMpc“^, the resulting confidence interval 
is found to be fag = 0.35 ± 0.17. 

While a significant improvement in precision is 
achieved by increasing the /Cmax value up to 0.3/iMpc“^, 
thereafter the gain is not as great as one might expect 
from the increased abundance of Fourier modes. One of 
the key limitations remains the degeneracy between the 
anisotropy parameter p and the velocity dispersion ay. 
Their joint likelihood is illustrated in Figure [m for two 
different values of /Cmax- Clearly if additional information 
were available to measure or predict the value of ay , sub¬ 


TABLE IV. The reduced values associated with each of 
the maximum likelihood values presented in Table |ml 


stantial improvements in the measurement of fag could 
be made. Another factor which limits the gains available 
from smaller scales is the shot noise. Its fractional im¬ 
portance is amplified by the reduction in the amplitude 
of the power, which becomes particularly apparent at the 
lowest threshold. 


VIII. DISCUSSION 

We have developed the clipping procedure proposed in 
Simpson et al [20] to enable its application to anisotropic 
fields, and applied this new analysis technique to the 
^ < 0.25 sample from the GAMA survey. A simple model 
based on the linear power spectrum is used to measure 
the normalised growth rate fag(z = 0.18) = 0.26 ± 0.10 
at /Cmax = 0.2 hMpc~^. Employing a higher order model 
from perturbation theory allows the use of higher wave 
numbers, while still not requiring numerical simulations 
for calibration. For this case we find fag{z = 0.18) = 
0.29 ± 0.10 when using /Cmax = 0.3/iMpc~^ and density 
fluctuations Sg < 4. These results alone are not in signifi¬ 
cant tension with expectations from the Planck data [43] 
within the context of a standard ACDM model. How¬ 
ever, they do add to a growing body of evidence that 
appears to prefer a lower amplitude of density perturba¬ 
tions at low redshifts. Such evidence includes weak grav¬ 
itational lensing [44], galaxy clusters [45], and a number 
of other measurements of redshift space distortions [46] . 
This trend is also visibly apparent in Figure but there 
are several possible explanations for this behaviour. One 
interpretation of this is the reduction in the quadrupole 
generated by nonlinear motions, relative to the Kaiser 
prediction, as illustrated in Figure 2 of [35]. However 
our result would be largely insensitive to this effect. An¬ 
other interpretation is the presence of nonlinear galaxy 
bias. An additional isotropic contribution to the power 
dilutes the strength of the anisotropic clustering signal. 
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FIG. 10. Constraints on fag from the power spectra of the galaxy field from GAMA. Results from the linear model are 
presented in the left hand panel, while those from the nonlinear model are illustrated in the central panel. The right hand 
panel uses the same model as the central panel but only includes intermediate wavevector angles (0.2 < /i < 0.8). Within each 
panel the squares, triangles and inverted triangles correspond to clipping thresholds (5o of 8, 5, and 4 respectively. The circles 
represent the maximum likelihood points when using the original field without clipping. 



FIG. 11. The joint likelihood contours at 68% and 95% when 
evaluating the GAMA power spectra with (5o = 4, using the 
linear model. The degeneracy with the velocity dispersion ay 
limits the gain in measuring the anisotropy parameter [3 when 
utilising higher wave numbers. The dashed contours show the 
degeneracy when fitting /cmax = 0.3/iMpc“^, while the solid 
contours represent /cmax = 0.5/iMpc“^. 


This effect can be seen in Figurewhere, before clipping 
is applied, the inferred value of fag is significantly lower 
than the correct value even when using k < 0.1 hMpc~^. 
The simplest forms of nonlinear bias are strongly sup¬ 
pressed by clipping, but others such as stochastic bias 
are likely to remain, and therefore merit further investi¬ 
gation. It is also important to note that most studies of 
redshift space distortions rely upon a prior on the range 


of possible background geometries based on results from 
WMAP. The portion of the error budget associated with 
the Alcock-Paczynski effect in each survey will therefore 
be highly correlated [ISIIITIIIH]. 

As was found to be the case in real space [9], the pre¬ 
ferred range of clipping thresholds are typically those 
that reduce the linear power by around 25 — 50%. For our 
galaxy field this corresponded to thresholds in the range 
4 < (^0 < 8- Higher thresholds lead to weaker clipping, 
which is less effective at removing the problematic con¬ 
tributions from nonlinear structure. Meanwhile stronger 
clipping from lower thresholds leads to reduced signal-to- 
noise, and may also induce a significant cross-correlation 
between the Gaussian and residual fields. 

At present the precision of our measurement of the 
growth of structure is limited by the degeneracy with 
the velocity dispersion ay. By applying a group finding 
algorithm to the galaxy catalogue it may be possible to 
reduce the influence of the Fingers of God. This would 
also improve the efficiency with which clipping removes 
peaks in the density field. 

Previous measurements of redshift space distortions ei¬ 
ther rely heavily on calibration from numerical simula¬ 
tions, or on more complex approaches to perturbation 
theory. Each of these rely on certain model-dependent 
assumptions, such as a linear bias model. Clipping is a 
complementary approach as it can identify whether the 
galaxy bias is showing signs of scale-dependence. The 
degree of covariance between parameter constraints ob¬ 
tained from the clipped analysis and a conventional non¬ 
linear analysis has yet to be quantified, however poten¬ 
tially these two approaches could be combined to yield 
considerable additional information within the same sur¬ 
vey volume. It may also be beneficial to perform a com- 
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bined analysis of power spectra from multiple clipping 
thresholds. 

Clipping may also be applicable to a number of other 
cosmological fields, which we shall consider in turn. 
Cosmic Microwave Background: 

In the early Universe, cosmological perturbations are 
know to be highly Gaussian. A recent analysis of the Cos¬ 
mic Microwave Background (CMB) from Planck severely 
limits the amplitude of local departures from Gaussian- 
ity /nl = 2.5 ±5.7 [49]. Applying a clipping transforma¬ 
tion to the CMB anisotropies is unlikely to be beneficial, 
since the Gaussian component is already highly domi¬ 
nant. However it may be of interest to identify whether 
features such as the lack of power on large angular scales, 
and the power asymmetry on the sky, remain intact after 
clipping, or are exacerbated. 

Lyman-a forest: 

At later times, the cosmological perturbations are again 
detectable in the absorption lines of quasars. This tech¬ 
nique has been used to detect the baryon acoustic oscil¬ 
lations in BOSS [50|. In this case the observed tracer 
already experiences a transformation similar to clipping, 
in that the highest density regions form damped Lya 
systems. The results of Section [TT| may be generalised to 
other local transformations. For example the transfor¬ 
mation relating the local density to the observed flux F 
is often approximated as 


g{S) = exp [-A(l + (5)“] , (38) 

which can be used in conjunction with <© to directly 
compute the flux correlation function in terms of the lin¬ 
ear power spectrum. 

Weak Gravitational Lensing: 

The shapes of high redshift galaxies are coherently dis¬ 
torted by the intervening matter perturbations. Cosmic 
shear offers the most direct insight into the dark matter 
distribution at lower redshifts, yet uncertainties in the 
small scale power spectrum limit the amount of cosmo¬ 
logical information that may be extracted. Gaussiani- 
sation of the convergence field have been proposed by 
several authors, such as Seo et al. [60], Joachimi et al. 
m, Yu et al [62] . However as highlighted in m, in the 
presence of shape noise the benefits of the transformation 
are minimal. This is due to the substantial reduction in 
the amplitude of the resulting power spectrum. Apply¬ 
ing clipping here with a suitably high threshold may be 
advantageous as it can suppress the strongest sources of 
nonlinearity while still preserving a high level of signal 
to noise. But the observable field appears in a projected 
two-dimensional form, due to the broad lensing kernel, so 
the identification and suppression of peaks is a less effi¬ 
cient procedure compared with the full three-dimensional 
data that can be acquired from the distribution of galax¬ 
ies. 

Galaxy Clustering: 

Cosmological information from the galaxy power spec¬ 
trum can be split into three categories: geometric infor- 
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FIG. 12. Comparison of measurements of fag as measured 
from anisotropic clustering in galaxy redshift data. These 
include 6dGFS [51]; WiggleZ [52]; BOSS [53] [54]; BOSS P(k) 
[55]; VVDS [56]; SDSS-II [57] [58]; VIPERS [59]; GAMA [^. 
The dotted and dashed lines represent theoretical predictions 
for a flat LCDM cosmology, as given by [Urn, erg] pairs from 
the maximum likelihood values from Planck @3] using only 
polarisation [0.286; 0.796], and the full dataset [0.3156; 0.831]. 


mation from the baryon acoustic oscillations; primordial 
information from the broader shape of the power spec¬ 
trum; and gravitational information from the degree of 
anisotropic clustering. 

Local density transformations such as clipping amelio¬ 
rate nonlinearities associated with high density regions. 
One form of nonlinearity for which this is not the case is 
that associated with the smoothing of the baryon acous¬ 
tic peak in the galaxy correlation function. Non-local 
transformations are more appropriate for this form of 
peak reconstruction, as demonstrated in [63U65] . where 
the signal can be largely restored by reversing the inferred 
large scale displacements. 

The large-scale shape of the galaxy power spectrum is 
sensitive to a variety of cosmological parameters such as 
the matter density, the spectral index, and the neutrino 
mass. The precision of these parameter measurements is 
limited by the uncertain nature of galaxy bias, which is 
expected to be linear on very large scales /c < 0.1 hMpc“^ 
but not on smaller scales where the vast majority of the 
information resides. Glipping can greatly assist in lin¬ 
earising the galaxy bias, thereby ensuring the clipped 
galaxy power spectrum bears a close resemblance to the 
clipped dark matter power spectrum. Interpreting the 
shape is less straightforward since stronger clipping leads 
to a change in the shape of the power spectrum, but this 
is fully specified by the transformation defined in (§. 
Overall, then, we see considerable scope for further ap¬ 
plications of the method presented here. 
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Appendix A: Spectral Distortion 

First we view how rapidly the series expansion given 
by ^ converges on the numerical solution. Taking the 
z = 127 snapshot as our fiducial Gaussian Random Field 
(GRF) we apply clipping and evaluate the true clipped 
power spectrum, Pc{k). We then generate an estimate 
of the clipped power, Pcjk), by truncating the series of 
^ at n = 0,2,4. Figure illustrates the ratio of the 
true and estimated power spectra in each case. In order 
to establish sub-percent precision in the estimated power 
spectrum, it is sufficient to stop at n = 2 provided the 
parameter ul > 0.5. Throughout this work we evaluate 
terms at n < 4. 



FIG. 13. The convergence of series given in to the nu¬ 
merical solution, as the number of terms in the summation is 
increased from n = 0,2,4. For this case the sample GRF was 
clipped at the 0.7(7 level, resulting in an approximately 40% 
loss of power on large scales. 


whose power spectrum recovers the standard form 


p{k,^l) = p%k){l + p^i^) 


2\2 


1 


1 + k^[pG^l2 


(B2) 


In Figure we can see that the fractional change in 
the angle-averaged power spectrum induced by clipping 
is slightly reduced when the velocity dispersion is in¬ 
troduced. The larger gradient in P{k) leads to a stronger 
contribution from higher order terms in However in 
all cases the dashed line of the model successfully repro¬ 
duces the behaviour of the data points. 

Next we evaluate the power spectrum from the mock 


galaxy fields at lower redshifts. Figure repeats the 
analysis of Figure but now with a smoothing length of 
I.95Mpc/h. The sets of parameter constraints derived 
from these two different smoothing lengths are found to 
be highly consistent with each other. 


Appendix B: Clipped Anisotropic Fields 

In this section we explore the consequences of Q by 
performing numerical tests on anisotropic fields. We take 
the 2 : = 127 Millennium-I density field in real space and 
impose a distortion along one axis consistent with the 
prescription of [26] : 


Appendix C: The GAMA Regions 


Figures 16 and 17 present constraints on /3 and fas 
from the three separate GAMA regions. In each case the 
results appear consistent between the three regions, and 
across the three thresholds within each region. 

The clipping statistics for each field are presented in 
Table El 






1 

yi + k‘^ix^al/2 ’ 


(Bl) 


Appendix D: Anisotropic Clustering 


where fi is the cosine of the angle between the wave vector 
and the line of sight, and the longitudinal amplification 
factor (3 = f /b. This leaves us with an anisotropic GRF 


Here we present constraints on the anisotropy parame¬ 
ter (3 from the clipped galaxy power spectrum. To do so 
we marginalise over the three model parameters (ul, cr^, 
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FIG. 14. The fractional change in t he a ngle-averaged power 
spectrum of the field described by ( |B2| , where we have set 
/3 = 0.5 and we explore (7 = 0 and cr = 5, each evaluated at 
two different clipping thresholds (removing 5% and 10% of the 
mass). The real-space density field is taken from the z — Vll 
snapshot of the Millennium-I simulation. Larger values of a 
can be seen to slightly reduce the effects of clipping on small 
scales. This arises from the contribution from terms in Q, 
which involve the self-convolution of the power spectrum. The 
horizontal dotted lines are the estimated clipped power when 
only using the first term from Q , while adding in the second 
leads to the dashed lines, which are in good agreement with 
the numerical results. 
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0.157 

15 
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0.024 

0.165 

9 

8 

0.005 

0.047 

12 

8 

0.008 

0.065 

15 

8 

0.008 

0.081 


TABLE V. The fraction of the field volume fv and the effec¬ 
tive fraction of the galaxies fm subject to clipping from each 
of the GAMA regions, for each of the three threshold values. 


bag) while fixing the shape of the linear power spectrum 
to the fiducial model. 

Figure illustrates the confidence intervals derived 
from the combination of the three fields. The maximum 
likelihood values for (3 are plotted as a function of the 
maximum wavenumber used to compare to the model. 
The behaviour closely reflects that of /ag, as seen in the 
central panel of Figure 
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FIG. 15. The same format as Figure]^ but we now use a finer grid to define the number density held, at 1.95Mpc//i instead of 
3.9Mpc//i. Since the total variance of the held is greater, slightly less power is lost after removing the same fraction of galaxies 
as before. 
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FIG. 16. Constraints on the anisotropy parameter ^ from the three individual GAMA regions, from left to right: G09, G12, 
G15, for the case of the unclipped galaxy field and three clipping thresholds (5o of 8, 5, and 4. 
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FIG. 17. Constraints on fag from the three individual GAMA regions, from left to right: G09, G12, G15, for the case of the 
unclipped galaxy field and three clipping thresholds (5o of 8, 5, and 4. 
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